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Abstract. We present a theoretical analysis of phase separation in the presence of a spatially periodic 
forcing of wavenumber q traveling with a velocity v. By an analytical and numerical study of a suitably 
generalized 2d-Cahn-Hilliard model we find as a function of the forcing amplitude and the velocity three 
different regimes of phase separation. For a sufficiently large forcing amplitude a spatially periodic phase 
separation of the forcing wavenumber takes place, which is dragged by the forcing with some phase delay. 
These locked solutions are only stable in a subrange of their existence and beyond their existence range 
the solutions are dragged irregularly during the initial transient period and otherwise rather regular. In 
the range of unstable locked solutions a coarsening dynamics similar to the unforced case takes place. For 
small and large values of the forcing wavenumber analytical approximations of the nonlinear solutions as 
well as for the range of existence and stability have been derived. 

PACS. 47.54.-r Pattern selection; pattern formation - 64.75.-g Phase equilibria - 05. 45. -a Nonlinear 
dynamics and chaos 



1 Introduction 

When a binary mixture at a critical composition is quenched 
^ from a homogeneous and usually high temperature phase 
. sufficiently below the coexistence curve then inhomoge- 
' neous concentration variations of a characteristic average 
, domain size occur which are known as spinodal decompo- 
' sition (SD) [I'2]. With the progress of time a coarsening 
process takes place, i.e. the average domain size increases 
continuously. The dynamics of phase separation under ho- 
mogeneous conditions has been intensively studied for nu- 
, merous systems during the recent decades |1I2| . In this 
> work we explore the effects of a traveling spatially peri- 
odic forcing on phase separation. 

For several kinetic mechanisms associated with first 
order phase transitions, such as nucleation and growth, 
thermal activation energy is required, but not for SD. The 
order parameter, that describes the phase separation in 
binary mixtures, is the local composition (concentration) 
and it obeys a local conservation law. Hence in phase sep- 
arating systems the order parameter is a conserved one in 
contrast to the large number of pattern forming systems, 
where one has to deal with non-conserved order parame- 
ters [3]. The phase separation is limited by diffusion- like 
processes and therefore in its late stage, where small in- 
homogeneities have developed into macroscopic domains, 
the decomposition dynamics becomes rather slow. 

Phase separation is also employed in producing techno- 
logically important functional elements in e.g. lithographic 
processes, biosensors and semiconductor devices ^4i5i6i7j . 
Therefore, an understanding of the mechanisms of the for- 
mation of regular structures is highly important for tech- 



nological applications and even the control of regular pat- 
terns is desired. 

In a large number of investigations the effects of exter- 
nal fields on phase separation have been examined, such 
as for instance the effects of shear flow jSj , time-dependent 
temperature variation [9] or local [TOj as well as spatially 
periodic heating [TT]. In the presence of spatially peri- 
odic forcing coarsening can be interrupted above a critical 
modulation amplitude and locked solutions can be found 

Especially for pattern forming systems with non-conserved 
order parameters external forcing has been recognized as 
a powerful method for investigations of the response be- 
havior of nonlinear patterns and of the inherently nonlin- 
ear mechanisms of selforganization. The effect of spatially 
periodic forcing on pattern formation has been investi- 
gated rather early in the context of thermal forcing [12] 
and electroconvection in nematic liquid crystals |13|14j . 
Further on, the response of stationary patterns with re- 
spect to static spatially periodic forcing in one or two 
spatial dimensions has been explored with respect to mod- 
ulated structures |13I14I15I16I17I18I19I20I2T] o r with re - 
spect to induced time-dependent phenomena |22|23|24] . 
Separately also the effects of temporal forcing have been 
investigated |25l26l27l28l29j . Recently a number of new 
phenomena have been found and explored by a combi- 
nation of spatial and temporal forcing in model systems 
30 31 32 33 34 35 36]. Here we present the first investiga- 
tion of the effects of pulled spatially periodic external forc- 
ing on a pattern forming model system with a conserved 
order parameter. 
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The paper is organized as follows. In Sec. 2 we present 
the model equation for the conserved order parameter 
field. In section 3 the one-dimensional periodic solutions 
and their existence range are described and their linear 
stability analysis in one and two dimensions is presented 
in section 4. The results of numerical simulations of the 
dynamics of phase separation under traveling spatially pe- 
riodic forcing in ID and 2D are given in Sec. 5. Some con- 
cluding remarks are added in Sec. 6. In the Appendix we 
describe the condition for the existence of spatially peri- 
odic solutions and analytical approximations for the solu- 
tions in the limit of small and large values of the forcing 
wavenumber. 



2 Model 

A model for phase separation in various systems is the 
Cahn-Hilliard equation for a conserved order parameter 
field ipivji) representing the local concentration of one of 
the components of a binary mixture. Here we focus on the 
effects of traveling spatially periodic forcing on phase sep- 
aration as described by the modified Cahn-Hilliard equa- 
tion 

dtijjir, t) = V^I-eiA + -0^ - + a cos[q{x - vt)]}. (1) 

The forcing term a cos[q{x — vt)] is caused by an interplay 
between a traveling temperature modulation and thermal 
diffusion (Soret effect). Its derivation has been described 
in detail for w = in Ref. [IT]. Such a traveling spatially 
periodic temperature modulation could be created for in- 
stance in optical grating experiments |37|38j with a light 
intensity /(x, t) ^ cos[q{x — vt)]. The control parameter e 
in Eq. ([1]) corresponds to a dimensionless distance to the 
critical temperature of the binary mixture. A transforma- 
tion of Eq. ([T]) in the frame comoving with the traveling 
forcing x ^ x — vt gives 

dti^{v, t) = [-£'0 + V'^ - V + a cos(ga;)] +vd^ip. (2) 



3 Periodic solutions and bifurcation diagrams 
in ID 

During the initial stage of an unforced phase separation, 
i.e. a — Q and v = Q, small perturbations of the homoge- 
neous solution = exhibit an exponential growth for 
e > and the growth rate cr = fc^(£ — fc^) shows a maxi- 
mum at a distinguished wavenumber km = \/e/2. Beyond 
the regime of exponential growth one observes the well 
known coarsening dynamics of phase separation. 

If phase separation is forced by a spatially periodic 
temperature modulation the coarsening dynamics is in- 
terrupted beyond a critical value of the forcing amplitude 
a and it is locked to the periodicity of the external forc- 
ing [11]. However, if this forcing is pulled by a velocity 
11 ^ 0, the spatially periodic solutions do exist only in a 
certain range of v depending on a. The conditions for the 



existence of spatially periodic solutions are formulated in 
Appendix I A . 1 1 and they are explicitely calculated for var- 
ious parameters in this section. 

The bifurcation diagrams of spatially periodic solu- 
tions of Eq. ^ are described for w = in subsection 13.21 
as well as for w ^ in 13.31 and the numerical method for 
calculating the solutions is explained in subsection 13. II 

3.1 Numerical solution 

The stationary 27r/(j'-periodic solutions of Eq. ([2]) can be 
represented by a truncated Fourier expansion 

N 

i^s{x)^ Cke'""" (3) 

k=-N 

with time-independent Fourier coefficients C_fe = CJl. Us- 
ing this ansatz one obtains after projection of Eq. ^ onto 

the j-th Fourier mode {J^^^'^ e~^'^^^ dx), the coupled sys- 
tem of nonlinear inhomogeneous equations for the corre- 
sponding expansion coefficients Cj 

Hjq)' - ijq)' + ^VJq]C, - ^{5,., + <5_i,,) 

-{jq)^Y.^'kCiCj-k-i^O, j^~N..N, (4) 

k,l 

which has been solved by Newton's method. N was ad- 
justed in order to keep a relative error smaller than 10^^ 
(e.g. iV = 81 for q = 0.5). 

3.2 Periodic solutions for v = 

For positive values of e spatially periodic solutions of Eq. ^ 
may be found with and without forcing. However, in the 
unforced case (a = 0) Eq. ^ has a ±V'-symmetry and 
therefore one has a pitchfork bifurcation from the trivial 
solution -0 = to finite amplitude periodic solutions as 
indicated by the dotted line in Fig. [1] For a = periodic 
solutions are unstable for any wavenumber q against coars- 
ening with the fastest growing mode as period doubling 

mm- 

For finite values of the modulation amplitude a the 
±?/'-symmetry is broken as indicated in Fig.[l]for a = 0.01 
by the dashed line and for a = 0.03 by the solid line. 
While we have in the unmodulated case a trivial solution 
■0 = and two finite solutions with identical amplitudes 
but of opposite sign, one finds in the forced case three 
periodic solutions, Ai, A2 and A3, of different amplitude 
as shown in Fig. [T] (b)-(d). The A^- and A2-solutions are 
in phase with the external modulation and the preferred 
Ai-solution is shifted by half a period. 

The amplitudes of the Ai- and the As-solution increase 
with increasing values of the control parameter e, but they 
also depend on q and a. For a given modulation amplitude 
a the amplitude of Ai decreases with increasing values of 



V. Weith, A. Krekhov, W. Zimmormann: Traveling spatially periodic forcing of phase separation 



3 




qx / (2jt) 



Fig. 1. Part (a) shows the bifurcation diagram for spatially 
periodic solutions at a forcing wavenumber q — 0.5 and mod- 
ulation amplitudes a = (dotted), a = 0.01 (dashed) and 
a = 0.03 (solid). Parts (b), (c) and (d) show the three so- 
lutions i^six) over one period corresponding to the different 
branches of the bifurcation diagram for a = 0.03 and e = 1. 

the modulation wavenumber q as indicated in Fig.[2l This 
trend holds also for the As-solution but the amplitude of 
the A2-solution increases with q. 

Keeping q fixed and increasing the modulation ampli- 
tude a the amplitude of the ^i-solution increases as can 
be seen in Fig. [TJ This also holds for the j42-solution. In 
contrast to this trend, the amplitude of the ^3-solution de- 
creases with increasing values of a (see dashed and solid 
lines in Fig.[T]). In case of small q (long wavelength modu- 
lation) the solutions Ai and A3 become steplike as shown 
in Fig. [5] (a) for Ai. This steplike shape also occurs, if the 
control parameter e becomes rather large compared to the 
modulation amplitude. The shape of the ^2-solution does 
not change when a or e are varied. 




0.5 1 0.5 1 0.5 1 

qx / (271) qx / (2n) qx / {2ti) 



Fig. 2. Ai-solution for different values of the modulation 
wavenumber: q = 0.25 (a), q — 0.5 (b) and q = 0.75 (c). The 
other parameters are e = 1, v = and a = 0.01. 



3.3 Periodic solutions for w 7^ 

The nonlinear solutions and the bifurcation diagram as 
given in Fig. [1] are only slightly changed by a finite travel- 
ing velocity. However, with increasing values of w a phase 



shift A(j) between the periodic forcing ~ cos{qx) and the 
solution ipsix) occurs. This is shown in Fig. [3] for three 
different values of v for both the Ai- and the Aa-solution. 
The ^i-solution, which is in antiphase with the external 
modulation, drags behind the forcing {A(p > 0) . This delay 
increases with increasing values of v (all other parameters 
fixed). The maximum phase shift that can be achieved is 
about — 7r/2 at a certain velocity v^x above which 

the solution does not exist. In contrast to this behavior 
the Aa-solution has a negative phase shift Acj) < with 
respect to the forcing as indicated by the lower part in 
Fig. [3l The maximum absolute value of the phase shift in 
this case is also tt/2. As indicated in Fig. |4](b) the critical 
phase shift A(j)m = tt/2 for the Ai-solution is reached with 
an increasing value of the velocity v at lower values of the 
control parameter e. 

For a fixed e a velocity Vex can be determined by the 
7r/2-criterion where both the Ai- and the As-solution van- 
ish to exist whereas the A2-solution still exists and propa- 
gates ahead of the external modulation with the absolute 
value of the phase shift always smaller than tt/2. However, 
as will be shown later, this solution is always unstable. 




0.5 1 0.5 1 0.5 1 

qx / {2n) qx / (2n) qx / (2n) 




0.5 1 0.5 1 0.5 1 



qx / (271) qx / (27i) qx / (27i) 

Fig. 3. j4i-solution (top) and yls-solution (bottom) for the 
velocity v^O (a,d), v = 0.008 (b,e) and v = 0.0142 (c,f). The 
other parameters are q = 0.5, a — 0.03 and e = 1. 

The amplitudes for the three solutions Ai [i = 1, 2, 3) 
are plotted as a function of the control parameter e for 
four different values of the propagation velocity v in the 
bifurcation diagram in Fig. HI As can be seen in this plot, 
the amplitudes for all four velocities are nearly identical. 

The shape of the solution ips becomes steplike for small 
values of q, as indicated in Fig. [2] (a). In the range q < 
0.3 the solution can be approximated rather well by a 
series of hyperbolic tangents as shown in Appendix IA.2I 
In the opposite case q > 0.5 one can achieve a rather good 
approximation of the periodic solution by restricting the 
Fourier ansatz in Eq. ^ to the leading modes as explained 
in Appendix IA.3I 
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Fig. 4. The bifurcation diagram for q — 0.5, a — 0.03 and 
finite values of the velocity v is shown in (a) and in (b) the 
phase shift A(j)/ {2tt) of the Ai-solution is presented for different 
velocities: v = 0.004 (solid), v = 0.006 (dashed), v = 0.008 
(dashed-dotted) and v = 0.0142 (dotted). 

4 Stability of periodic solutions in ID and 2D 

The stationary periodic solutions V's(a;) of Eq. ([2]) stud- 
ied so far are stable only in a subrange of parameters 
where they exist. In this section we investigate the sta- 
bility of periodic solutions ips{x) of Eq. ([2]) with respect 
to small perturbations. The resulting linear partial differ- 
ential equation for the perturbation V'l 

dt^x, t) = dl {-e + 3^-,' - dl) + vd.iJi (5) 

contains a spatially periodic coefficient depending on ips- 
Accordingly the solutions of Eq. ([5]) are due to the Floquet 
theory of the form 

V'i(x,t) =e^* e^^" (6) 

with the Floquet parameter s and a 27r/(7-periodic function 
(ppix). This periodic function (t>F{x) can be represented by 
a Fourier expansion 

N 

<Pf{x)= ^ne^'™- (7) 

n=-N 



Taking into account Eq. ([3|) for ^pg the linear partial dif- 
ferential equation (O is transformed into an eigenvalue 
problem 

aDi = [e{ql + s)2 - {ql + s)* + iv{ql + s)] Di 

-3((7/ + s)2^CfeC^A-fc-m, l^-N..N, (8) 
k,m 

with the coefficients Ck determined by Eq. (|4]). 

The growth rate a for the solutions Ai , A2 and is 
shown in Fig. [5] as function of the Floquet parameter s 
for three different values of the velocity v and fixed val- 
ues of the forcing amplitude a, the forcing wavenumber q 
and the control parameter e. We found that for the so- 
lutions A2 and A3 the growth rate a is positive for all 
combinations of the parameters e, q, a and v. Accord- 
ingly A2 and A3 are always unstable. The growth rate of 
the perturbation of the ^i-solution is negative for certain 
parameter ranges, meaning that this is the only solution 
that can be stabilized. For fixed e and q the modulation 
amplitude a has to exceed a certain value as{q) to stabi- 
lize the Ai-solution. If the traveling velocity v is smaller 
than a critical one Vc{s,a,q) the solution remains stable. 
The critical velocity Uc(£, a, q) is given by the solid line in 
Fig.mfor a certain set of parameters and for v > Wc(e, a, q) 
the periodic solutions are linearly unstable. The onset of 
instability occurs for small values of the Floquet exponent 
s 0, i.e. belongs to a long-wave perturbation. In Fig. [S] 
also the boundary of the existence range of periodic so- 
lutions Vex{s, a, q) is shown as obtained approximately by 
Eq. (|27|) (dashed line) and by a full numerical simulation 
(dotted line). Since the stability boundary (solid line) al- 
ways lies below Vf,x the periodic solution always becomes 
unstable before the existence range is reached. 
Extending the analysis to two spatial dimensions the sta- 
bility of a two-dimensional stripe pattern that is periodic 
in the x-direction can be investigated using ip{x,y,t) = 
ips(x) + ij^2ix,y,t) with 1^2! << IV'sl- The procedure is 
the same as in one dimension whereas the perturbation 
ip2{x,y,t) 

Mx,y,t)^e-' e'^'-+Py^ Mx), (9) 

depends now on two Floquet parameters s (x-direction) 
and p (y-direction). Linearizing Eq. ^ and using expan- 
sion ([7]) one obtains the eigenvalue problem 

aDi = [e{ql + s)2 + ep2 _ _^ ^^4 

~2{ql + s fp^ - / + iv{ql + s)]Di 

-?,[{ql + sf + /] C'kC'raDi^k-m , I = -N..N. 

The linear stability analysis shows that additional transver- 
sal degrees of freedom do not influence the stability dia- 
gram obtained for the one-dimensional case. 

5 Dynamics of the phase separation 

To study the temporal evolution of the phase separation 
numerical simulations of Eq. ([1]) for one and two spatial 



V. Weith, A. Krekhov, W. Zimmormann: Traveling spatially periodic forcing of phase separation 



5 




Fig. 6. Above the solid line the spatially periodic solution is 
unstable. The dashed line marks the existence boundary above 
which the spatially periodic solution does not exist due to the 
criterion given by Eq. (|27p . The dotted curve marks the ex- 
istence boundary obtained numerically from Eq. Q. The pa- 
rameters are q — 0.5 and e = 1. 
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Fig. 5. The growth rate a of small perturbations with respect 
to the stationary periodic solution is shown as a function of the 
Floquet parameter s for the Ai (a), the A2 (b) and the ^3- 
solution (c) for the three propagation velocities v — (dotted), 
V = 0.008 (dashed) and v = 0.014 (solid). Further parameters 
are q = 0.5, a = 0.03 and e = 1. In case of the yl2-solution all 
three lines are indistinguishable. 

dimensions have been performed. We have used a central 
finite difference approximation of the spatial derivatives 
and an Euler integration of the resulting ordinary differen- 
tial equations in time. Periodic boundary conditions have 
been used and we start initially with the homogeneous 
state 1^ = with small superimposed noise. 



by the structure factor at the modulation wavcnumbcr q 

S{q,t) = \M<lM , ^P{q,t) = J dxe''^-i>ix,t). (11) 

In the early stages of phase separation the amplitude of 
S{q,t) grows faster for w ^ and even exceeds its final 
value for some time as can be seen in Fig. [T] The tempo- 




Fig. 7. The dynamics of the maximum of the structure factor 
(relating to its final value) for v — (solid line) and v — 
0.0266 < Vc (dashed line). The other parameters are e = 1, 
q = 0.5 and a = 0.1. 



5.1 ID Simulations 

In one dimension we have discretized a system of length 
L = 24:TT / q hy N = 200 steps with a step size 6x = L/N 
and a time step At — 10^^. The noise amplitude for the 
initial condition — was adjusted at 10^^. The evolu- 
tion of the inhomogeneous pattern has been characterized 



ral evolution of ipix, t) is shown in Fig. [8] for three veloc- 
ity regimes. In part (a) of Fig.[8]the velocity is sufRciently 
small and belongs to the range where the pattern is locked 
to the propagating external modulation. In part (b) the 
velocity is chosen in the range Vc < v = 0.0185 < Vex 
where the solution locked to the traveling external modu- 
lation still exists, but where it is linearly unstable. In this 
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Fig. 8. Temporal evolution of phase separation is shown for 
V = 0.015 < Vc (a), Vc < V = 0.0185 < Vex (b) and v = 
0.02 > Vex (c,d). The other parameters are e — 1, a — 0.04 and 
g = 0.5. 



parameter range the solution is locked during the initial 
period of phase separation before coarsening takes over. 



In parts (c) and (d) the temporal evolution of phase sep- 
aration is shown for a velocity v = 0.02 > v^^ where the 
locked solution does not exist anymore. At this velocity an 
interesting pinning-depinning behaviour can be observed 
during the initial stage of phase separation. One still has a 
travelling periodic solution with the same wavenumber as 
the forcing but with a velocity smaller than the velocity of 
the forcing. Due to the velocity mismatch the phase shift 
between the solution and the forcing is slowly increased 
before it reaches about half of the forcing period. From 
that moment the periodic solution practically stops mov- 
ing (pinning) until the forcing shifts over the next half of 
the period. After that the solution starts moving again 
(depinning) and the process repeats itself a few times. 
Later on the agreement of the wavelengths of the solution 
and the forcing cannot be maintained anymore and the 
coarsening takes place. 



5.2 2D Simulations 

In two spatial dimensions the same system size as in 5.1 
has been used in both directions, = Ly = 24,71 /q with 
Nx = Ny = 256 discretization points. Again we chose lO""' 
as noise amplitude for the initial homogeneous state. 
The numerical simulations in two dimensions are in agree- 
ment with the results of the linear stability analysis of Sec. 
4. 




Fig. 9. Snapshot of the phase separation in 2d for q — 0.5, 
V = 0.0188 1.07 ■ Vc), a = 0.04, e = 1 at t = 5 x 10'' (a) 
and for q = 0.75, v = 0.266 (r^ 1.32 • Vc), a = 0.3, e = 1 at 
t = 5x 10^ (b). The curves at the bottom indicate the external 
forcing. 



For a > ac and v < Vc one obtains stable stripe patterns. 
In the unstable range v > Vc the locking of the solution 
gets lost and a coarsening sets in as shown in Fig. [9] (a) 
for V « 1.07 • Vc slightly beyond the critical velocity and 
in Fig.[9](b) for v = 0.266 (« 1.32 -Vc). For v > Vc the lin- 
ear stability analysis of the locked solutions in two spatial 
dimensions shows that transversal perturbations being pe- 
riodic in the y-direction exhibit also a positive growth rate 
(t(s,p). However, for small values of the forcing wavenum- 
ber q, the modulational wavenumber p is very small and 
the corresponding length 2tt/p may be larger than the 
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system size in our simulations. In this case the coarsen- 
ing is quasi-one-dimensional along the x-direction as can 
be seen for instance in Fig. [9] (a). For larger values of q 
also the length scale 2Tr /pc, with pc at the maximum of 
the growth rate a{s,p), becomes smaller than the system 
length Ly. In this parameter range the coarsening shows 
a two-dimensional behaviour as shown in Fig. [9] (b) . 



6 Summary and conclusions 

Effects of a spatially periodic forcing, which travels with 
velocity v, on phase separation in binary mixtures have 
been investigated in terms of a generalized Cahn-Hilliard 
model. For a sufficiently large forcing amplitude we identi- 
fied three characteristic velocity regimes which are distin- 
guished by the different evolving pattern and their dynam- 
ics. For a vanishing velocity v and the forcing amplitude 
a larger than a critical value Oc the phase separation can 
be locked to the external modulation [TT] . 

For a finite velocity v a condition is derived in Ap- 
pendix lA.ll which determines the parameter range where 
the spatially periodic solutions of the same wavelength 
as the external forcing exist. This spatially periodic so- 
lution is locked with respect to the external forcing in a 
finite parameter range beyond a > Qq and it is dragged 
by the forcing with some phase delay. This phase delay 
increases with the velocity v up to some critical value Vex ■ 
Beyond Vex the locked spatially periodic solution does not 
exist anymore as described in Section [33] and as shown by 
Fig.Sl 

In the limit of small values of the modulation wavenum- 
ber q the forced solutions resemble a periodic array of 
kinks and may be approximated by a series of hyperbolic 
tangents as given by Eq. (fT6| . For this approximate so- 
lution we derived in Appendix IA.2I with the help of the 
existence condition given in Appendix lA.ll an analytical 
expression for the boundary of the parameter range of ex- 
istence of the locked solutions. Also in the range of larger 
values of q an analytical approximation of the nonlinear 
solution as well as its boundary of the existence range has 
been determined in Appendix I A. 3 1 

The locked and dragged spatially periodic solutions 
are linearly stable only in a subdomain of their range of 
existence as described in Section S) The critical velocity 
Vc{e, a, g), beyond which the locked solutions become un- 
stable and a depinning of the locked solutions takes place, 
is given as function of a in Fig. For large values of q 
the numerical results of the stability analysis have been 
compared with results of an analytical approximation in 
Appendix EH cf. Fig. [Til 

In the parameter range of unstable locked solutions 
one observes in the two phase region a coarsening dy- 
namics rather similar to the coarsening in the unforced 
case. However, during an initial transient period the non- 
linear evolution of the solutions is different in the param- 
eter range where the locked solutions still exist from that 
where they do not exist anymore. In the former case they 
are just locked during a transient period to the traveling 
forcing and in the latter case irregular phase jumps occur. 



In the range, where the locked solutions are either un- 
stable or do not exist, we observe for finite values of the 
traveling velocity v a faster coarsening dynamics than in 
the limit v = 0. 

The analysis presented here shows that travelling spa- 
tially periodic forcing can serve as an effective tool for the 
control of phase separation. The finite velocity of the forc- 
ing leads to a faster dynamics of the initial stage of phase 
separation compared to the case of a stationary spatially 
periodic forcing. The results should be relevant to thin 
polymer blends with a moderate Soret coefficient and a 
weak periodic temperature modulation created, e.g., by 
means of optical grating techniques. 



Acknowledgments 

We thank C. Feller, W. Kohler and W. Pesch for enlight- 
ening discussions. Financial support by DFG through SFB 
481 is gratefully acknowledged. 



A Stationary periodic solutions 

A.l Existence of spatially periodic solutions 

The condition of existence of stationary periodic solutions 
ips{x) of Eq. ([2]) with the wavenumber q of the external 
forcing is derived in this section, ijjs (x) is a solution of the 
following equation 

= dl[-e^s +'>Ps- dl^s + acosiqx)] +vdx^s- (12) 

Assuming 4's{x) being a periodic function with a vanish- 
ing mean value /g"/'i^,(x')d:E' = one obtains after two 
integrations of Eq. (|T2|) 



dx <! -e^V'f + - ^{dxips)'^ } + acos{qx){dxtps) 



+v (dxips) J 'ips{x')dx' = C {dxi^s)- 





(13) 



Using again the periodicity of V's, an integration of 
with respect to the interval [0, 27r/(7] gives 



acos{qx){dxi}s) + v(dxi^s) / ipsix')dx' 



2-K/q 

0= J dx 

_ _ 

(14) 

After an integration by parts of the constraint of ex- 
istence for stationary periodic solutions ips{x) of Eq. ([2]) 
follows: 
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A. 2 Analytical approximation in the limit of small q 

For small values of the wavenumber q the solutions for 
the Ai-branch become steplike periodic functions as indi- 
cated in Fig. [2] (a) . Kink-type periodic functions may be 
approximated by a series of hyperbolic tangents ^SQj 



JV 



ipsix) ^ ^/i-^{-lTt&nh [M{x)] 



n=0 



where 



M{x) 



n+- -+7 
2 q 



(16) 



(17) 



describing a periodic array of N single interfaces with a 
spacing Tr/q (size of the system L = N ■ n/q). 7 depends 
on q, e, a and v. Inserting the ansatz p6p into Eq. (jlSp 
one obtains 



in terms of the parameters of the external forcing. With 
7 as determined above the representation in Eq. (|16p is 
a good approximation of the solution as it deviates less 
than 10% from the full numerical solution for q < 0.3. 
The deviation decreases with decreasing q. 

In the range of the coexistence of the three branches 
Ai, A2 and A3 the shift 7 between the forcing and the 
spatially periodic solution increases with an increasing ve- 
locity V up to the numerically determined existence limit 
of the yli-solution at approximately jq — A(j)„i = 7r/2, as 
described in Sec. 13.31 With this condition jq = tt/2 we 
may also calculate the velocity 



2aq 



(23) 



below which the periodic kink-type solution exists. This 
analytical expression for Vex is in a good agreement with 
the numerically determined existence boundary. 



27r/g 



/ sin(q2;) • ^/i ^(-1)" tanh [M(x)] da; 



27r/g 



V 

aq 



N 



Vi^(-l)"tanh[M(a 



Tl = 



dx (18) 



for the determination of 7, which may be reformulated in 
the following form 



/ sm{qx) ■'^{-iyH'Anh[M{x)]dx 



w^e j 27r 
aq\q 



/ ^sech^ [M(x)]da; I. (19) 

n n=0 J 



For g <C 1 the first integral in this equation may be ap- 
proximated by the expression 

/ sin(ga;) • ^(-1)" tanh [M(x)] dx « - sin(79) 



and the second integral by 

l-Kjq 



/ ^sech2[M(x)]da;«4W-. 



(20) 



(21) 



Both integrals together with Eq. give an explicit ex- 
pression for the constant 7 



1 

7 = — arcsm 
9 



V\/£ / TT 



a \2q \ e 



(22) 



A. 3 Analytical approximation in the limit of large q 

For large q the nonlinear solution can be described by a 
leading order Fourier expansion: ipg cx cos{qx). However, 
if the forcing is propagating with a velocity v a phase 
shift between the forcing and the periodic solution ips(x) 
occurs and this phase shift can be taken into account by 
the following ansatz: 



tpsix) = Ci cos{qx) + C2V sin(ga;). 



(24) 



This corresponds to a restriction of the series given by 
Eq. jS]) to TV = 1. Inserting the ansatz ([24|) in Eq. (fT2l) 
and collecting the linear independent contributions pro- 
portional to sn\{qx) and cos{qx) one obtains the two cou- 
pled equations for the coefficients Ci and C2 



s-q^-\cl-^-^Cl]c,+-C2-a = Q, (25) 



e~q'~\cl-\ci\c,--C,^Q. (26) 



In the range q > 0.5 a numerical solution of these two 
coupled equations provides solutions which deviate less 
than 30% from the full numerical solution. 

The deviation decreases further with increasing q. A 
comparison between the solution ((24|) and the full numer- 
ical solution is shown in Fig. \W\ for q = 0.5. The devia- 
tions increase at larger values of e and close to the boarder 
where the stationary solutions vanish to exist. The exis- 
tence range of the ^i-solution again may be determined 
by assuming a maximum phase shift of 7r/2 between the 
external forcing and the periodic solution. In this case 
the cos-contribution in Eq. (|24p vanishes, that means one 
should have Ci = 0. Then the velocity Vex can be calcu- 
lated from the existence of the solution C2 7^ of Eqs. ([25|) 
and jag): 

- ^ "'^ :. (27) 
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1.5 




0.5 1 1.5 2 



E 

Fig. 10. The bifurcation diagram obtained from the full nu- 
merical solution (solid) and the approximation (|24p (dashed) 
for the parameters q = 0.5, a = 0.03 and v = 0.01424. 



>" 0.15- 




Ql , ^ , , , 1 

0.34 0.36 0.38 0.4 0.42 0.44 0.46 
a 

Fig. 11. Critical velocity from the numerical stability analysis 
(solid) and the analytical approximation (|25[) . (|26|l and (|29p 
(dashed) for q = 0.9 and e = 1. 



The comparison with the numerical boundary is shown in 
Fig.l 

In the limit of large q it is also possible to obtain an an- 
alytical approximation for the linear stability of the locked 
periodic solutions. Therefore we use Eq. ([6]) as an ansatz 
for the perturbation ipilx^t). Now we approximate (t>F{x) 
by the three dominant Fourier-modes 

= i^„2e-2*«" +Do + 026^'"'=. (28) 

Inserting ([M)) and (pS)) into Eq. ([5]) we gain an eigenvalue 
problem. Then we expand the growth rate a in terms of 
the Floquet parameter s and use the curvature for s ^ 
as a criterion for the stability boundary, this leads to 

= 27g2c|?;6 - ISOq'^CfC^v^ - 256q'^ + 128q^ - Av^ 
-16(7^ - 384q'^Cf + 72q^Cf - 90q'^Cf - iSAq^C^v^ 
-90q'^C^v^ + 12q^Clv^ + iMq^Clv^ + 2\Qq'^Ctv'^ 
+?,8'iq^Cl + 8lq^Ctv'^Cl + Slq^CfCiv^ 
+432q^Civ^Cf + 6Civ^ + QClv^ -t- 216q^C^ 
+21q^Cl (29) 

as a condition for linear stability. Together with Eqs. 
(PS)) this gives a system of three coupled equations that 
can be solved numerically. In Fig. [TT] this approximation 
is compared with the results of the full numerical stability 
analysis. 
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